Relativistic polarization analysis of Rayleigh scattering by atomic hydrogen 
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A relativistic analysis of the polarization properties of light elastically scattered by atomic hydro- 
gen is performed, based on the Dirac equation and second order perturbation theory. The relativistic 
atomic states used for the calculations are obtained by making use of the finite basis set method 
and expressed in terms of B splines and B polynomials. We introduce two experimental scenarios 
in which the light is circularly and linearly polarized, respectively. For each of these scenarios, the 
polarization-dependent angular distribution and the degrees of circular and linear polarization of 
the scattered light are investigated as a function of scattering angle and photon energy. Analytical 
expressions are derived for the polarization-dependent angular distribution which can be used for 
scattering by both hydrogenic as well as many-electron systems. Detailed computations are per- 
formed for Rayleigh scattering by atomic hydrogen within the incident photon energy range 0.5 to 
10 keV. Particular attention is paid to the effects that arise from higher (nondipole) terms in the 
expansion of the electron-photon interaction. 

PACS numbers: 32.80.Wr, 32.90.+a 



I. INTRODUCTION 



Polarization is one of the main characteristics of light 
which can be employed in order to investigate the prop- 
erties of matter. In atomic physics, especially, the po- 
larization properties of light have been studied for vari- 
ous processes, such as atomic and ionic photoionization 
1 1 , hyperfine-quenched transitions I2j , two-photon decay 
3 , H[ , atomic- field bremsstrahlung [5[ , radiative electron 
capture [gj and clastic scattering of light. Elastic scat- 
tering (so called Rayleigh scattering) by atoms and ions 
0, [H has applications in astronomy, shielding, medical 
diagnostics and also is used extensively to obtain infor- 
mation about the structural properties of materials and 
complex molecules ■ 

Owing to the develepment of x-ray polarization sensi- 
tive detectors @, tunable polarization free-electron 
lasers [l6[ as well as synchrotron radiation sources [l7j |. 
an increasing demand for accurate theoretical prediction 
on polarization-dependent atomic phenomena has been 
pointed out in the literature. For instance, elastic scat- 
tering experiments with a nearly 100% polarized hard X- 
ray beam have been recently performed in ESRF (Greno- 
ble, France), while further experiments are planned to 
be realized at the DESY synchrotron facility (Hamburg, 
Germany) in the near future [l8|. Moreover, the novel 
linear polarimetry technique based on Rayleigh scatter- 
ing was used for the first time [l5| and recently applied 
as a complementary techninque to Compton scattering 
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During the past decades, the total cross section for 
Rayleigh scattering has been widely investigated within 
the relativistic as well as nonrclativistic frameworks 
[25j . In contrast, only a few studies have been done on the 
angular and polarization properties of scattered photons. 
For example, Roy et al. calculated linear polarization ef- 
fects in the elastic scattering of x-rays and 7-rays from 
targets with Z ranging from 13 to 92, by using numeri- 
cally obtained Rayleigh scattering amplitudes [26[ . More 
recently, Manakov et al. analyzed photon-polarization 
effects in two photon bound-bound atomic transitions, 
including relativistic and retardation effects p7| . In par- 
ticular, they numerically investigated circular dichroism 
effects in scattering of hard photons by hydrogenic as 
well as many-electron systems 0, [29j|. However, to the 
best of our knowledge, a fully relativistic treatment of the 
polarization-dependent angular distribution in Rayleigh 
scattering by atomic hydrogen has not been performed 
yet. 

In the present work, we investigate the polarization 
properties in Rayleigh scattering by an unpolarized hy- 
drogen atom, within the fully relativistic framework of 
the Dirac equation. We introduce two experimental sce- 
narios in which the light is circularly and linearly po- 
larized, respectively. For each of these scenarios, we in- 
vestigate the polarization-dependent angular distribution 
(PDAD) and the degrees of circular and linear polar- 
ization of the scattered light. We derive an analytical 
expression for the PDAD, which is valid for scattering 
by hydrogenic as well as many-electron systems. With 
the aid of the Wigner-Racah algebra, we write such ex- 
pression in terms of angular parts and reduced matrix 
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elements, where the latter are independent of the scat- 
tering geometry. The numerical evaluation of the re- 
duced matrix elements is relativistically carried out for 
the Rayleigh scattering by atomic hydrogen, through the 
use of finite basis sets for the Dirac equation constructed 
from B splines and B polynomials. The photon energy 
range we investigate is 0.5 to 10 keV. It has been re- 
cently showed by us that, within this energy range, the 
finite basis set approach to Rayleigh scattering gives good 
agreement with NIST (National Institute of Standards 
and Technology) data values and other calculations [3(| ■ 
This article is structured as follows. In Sec. [TT1 we in- 
troduce the geometry and the notation used. In Section 
IIII1 we recall the general polarization-dependent transi- 
tion amplitude for Rayleigh scattering and evaluate it 
separately for the circular and linear polarization scenar- 
ios. In Section IIV1 we describe the numerical method 
used for carrying out the calculations, while, in Sec. [Vj 
we present the results for the PDAD and the degrees of 
circular and linear polarization of the light scattered by 
atomic hydrogen. Finally, a short summary is given in 
Sec. EH 

SI units are used throughout the article. 



axis along the direction of the incident photon (fei). As 
we will see, such a choice of quantization axis simplifies 
the multipole expansion of the electron-photon interac- 
tion operator. The scattered photon propagates along 
the direction k 2 at angle 9 with respect to the z axis. 
The scattering plane (xz) is defined by the incident and 
scattered photon's directions. Incident (scattered) pho- 
ton has energy S 7l , 2 . = fiw 1 ( 2 ), propagation vector fci(2) 
and polarization unit vector £i(2) , where H is the reduced 
Planck constant. 

The photon polarizations we consider are linear (i = 
e x ) and circular (e = e^). Neither linear combinations 
of these polarizations nor mixed polarization states are 
taken into account. In particular, we do not consider 
clliptically polarized light. From an experimental point of 
view, linear polarization is probably the most interesting 
type for the considered photon energy range [fjj]. The 
polarization angle < \ < n and the helicity A = ±1 are 
the variables used to parametrize the linear and circular 
polarization states, respectively. In this notation, A = +1 
describes right-handed and A = — 1 left-handed polarized 
photons, respectively. The definitions of the polarization 
angles are displayed also in Fig.[T]if incident and scattered 
photons are taken to be linearly polarized. 



II. ATOMIC SYSTEM AND GEOMETRY 

Let us start by introducing the atomic system and the 
geometry under which the distribution of the Raylcigh- 
scattered photons is investigated. We consider hydrogen 
atom in the ground state which is irradiated by light, 
as displayed in Fig. [TJ We adopt the quantization (z) 
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FIG. 1: (Color online) The adopted geometry for the scatter- 
ing process is displayed. The polar angle 9 uniquely defines 
the direction of the scattered photon in the xz plane (scat- 
tering plane). The angle xi (x^) parametrizes the linear po- 
larization of the incident (scattered) photon. The hydrogen 
atom is placed at the origin of the coordinate axes xyz. 



III. THEORY 

A. Transition amplitude and observables 

The second-order transition amplitude for Rayleigh 
(and Raman) scattering is given by [30l - l32T | 



Mif{ei,i 2 ) 



(/|ftt(fc 2 ,e 2 )H Hft(fci,6i)|i) 



t 



(/Iftffc!,^)!^) (v\rt(k 2 ,i 2 )\i) 



, (1) 



where u V i = (E v — Ei)/h is the transition frequency be- 
tween states \v) and \i). Here, the transition operator 
^•(^1(2)1 ^1(2)) describes the relativistic interaction be- 
tween the bound electron and the incident (scattered) 
photon. In the Coulomb gauge, the explicit expression of 
this transition operator is 



lZ(k, e) — a. ■ e e 



ih-r 



(2) 



where a. is the vector of Dirac matrices. 

The summation over the intermediate states showed 
in Eq. ([T]) runs over the complete one-particle spectrum 
| v) , including a summation over the discrete part of the 
spectrum as well as the integration over the positive and 
negative energy continua. Performing this summation is 
perhaps the most difficult part in determining the transi- 
tion amplitude and we shall postpone further discussion 
on it to Sec. [TV] 

The initial \i) and final |/) states of atomic hydrogen 
have well-defined angular momentum j, angular momen- 
tum projection rrij and parity (—1)', where I is the orbital 
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angular momentum of the larger component of the Dirac 
spinor. In the following, we denote them respectively as 
jijiriji) and \/3f,jf,mj.}, where /3 is a collective la- 
bel used to denote all the additional quantum numbers 
needed to specify the atomic states but for j and rrij. 
For hydrogenic ions, specifically, (3 refers to the principal 
quantum number n and the parity quantum number 
Owing to the conservation of energy, moreover, the ener- 
gies i? 7l 2 and Ef^ are simply related by the equation 



E f - Ei 



E~ 



Erv 



(3) 



Since Rayleigh scattering is an elastic process, the initial 
and final states coincide, |i) = |/), and, thus, Eq. (|3]) 
simplifies to E 7l = E l2 = E~ r It is furthermore assumed 
that both initial (m^ ) and final (rrij, ) polarizations of 
atomic states remain unobserved, as typical for most ex- 
periments. 

In this work, we shall separately consider the two ex- 
perimental scenarios corresponding to circularly and lin- 
early polarized light. In the circular polarization sce- 
nario, the incident light is circularly polarized or un- 
polarized and the polarization of the scattered light is 
measured in the circular base. In the linear polarization 
scenario, the incident light is linearly polarized or un- 
polarized and the polarization of the scattered light is 
measured in the linear base. We conduct our analysis by 
investigating the PDAD (i.e., the polarization-dependent 
angular distribution) and the degrees of circular and lin- 
ear polarization of the scattered light. 

For the circular polarization scenario, the PDAD can 
be written in terms of the scattering amplitude (HJ as 

M 



da c 
dO 



(Ai, A 2 , E n , i 



a 2 c 2 
(2j< + 1) 



|-^ c (Ai,A 2 ) 
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(4) 

if the polarization of incident light is known. In the above 
equation, a is the electromagnetic coupling constant, c is 
the speed of light in vacuum and for simplicity we defined 



M c (\ 1 ,X 2 )=M if (e c Xl ,et 2 ) . 



(5) 



If the incident light is unpolarized, the PDAD for the cir- 
cular polarization scenario is obtained by taking (0]) and 
averaging over the (two independent) circular polariza- 
tions of the incident light [33| : 



da c 
dO 



(A 2 ,£. 



7) ' 



^E^( A i> A ^7,0)- (6) 



The PDADs for the linear polarization scenario, 
da 1 da 1 

-^(xi,X2,E 1: 9) and — (x2,E 1 ,9), are simply ob- 
tained from Q and ((6|) respectively, with the replace- 
ments c — > I and Ai,2 — > Xi,2- 

The degrees of circular (Pc) and linear (Pl) polariza- 
tion of the scattered light, if incoming light is polarized, 



are defined as @, Hi| 

Pc = Pz 
Pl 



(P2) 



(7) 



where 
Pi = 



da' (xu 0°, E 7 , 9) - dg'Qq, 90°, £ 7 , 1 
da l {xi^°,E 1 ,e)+da l {xi^°,E 1 ^ 



= da'(xi, 45°, E 7 ,9)- da l ( X i, 135°, E 7 , 9) (g) 



Pa 



da 1 ( X iA5°,E 1 ,9) + da l (xi, 135° , E 
dcr c (Ai,+l,^ 7 ,6») - do- c (Ai,-l,£ 7 ,l 



7 ' 1 



do- c (Ai, +1, E y , 9) + da c (Ai, -1, E y , 9) ' 

Pi, 2. 3 are called the first, second and third Stokes pa- 
rameter, respectively. We have denoted, above and in 
the following, da/dfl by da for brevity. If the incom- 
ing light is unpolarized, the degrees of circular and linear 
polarization of the scattered light will be denoted by Pc 
and Pj,, respectively Their definitions arc given by 



Pc = h , 

Pl = J(Pi¥ + (P 2 y 



(9) 



where 



- d(7^ (0° , E*y , 9) dfj' (90° , Ery , 9) 

1 ~~ da l (0°,E 7 ,9) +da l (90°,E 7 ,9) ' 



Pi = 



Pa 



dg'(45°,£ 7 ,6Q -dg'(135°,£ 7 ,6>) 
da l (Ab°,E 1 ,9)+da l (l?>5 ,E 1 ,6) 

da c (+l,E 7 ,9) -dCT c (-l,£ 7 ,6>) 



(10) 



da c (+l,E 1 



da c (-l,^ 7 ,( 



Prior to showing the results for the above defined 
PDADs and degrees of linear and circular polarization of 
the scattered light, in the following subsections we will 
further evaluate the scattering amplitude ([TJ for both 
circular and linear polarization scenarios. We will show 
that the use of Wigner-Racah algebra technique allows 
for significant simplifications. 



B. Evaluation of the transition amplitude for 
circular polarization scenario 

We shall here evaluate the amplitude ([T]) for the case 
in which the incident light is circularly polarized or un- 
polarized and the polarization of the scattered light is 
measured in the circular base. To this end, we expand 
the vector plane wave e x e lkr in terms of spherical tensors 
with well defined angular momentum properties j35| : 



pc ik-r 



^E E E i L [L] 1/2 {i\Ya? LM {k,r) 

L=\ M=-Lp=0,l 

xD L MX ^ k ,9 k ,0) , (11) 
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where [L 1} L 2 , L n ] = (2L, + l)(2i 2 + l)...(2L n + 1), 
and the spherical tensor a p LM (k,r) refers to the mag- 
netic (p = 0) and electric (p = 1) multipolcs. Each term 
a p LM (k, r) has angular momentum L, angular momentum 
projection M and parity (— 1) L+1+P . 

As seen from Eq. (fTT|) , the angular dependence of the 
vector plane wave results from the Wigncr (rotation) ma- 
trices. The Wigner matrices transform each multipolc 
field, with original quantization axis along the photon 
propagation direction, into the field with quantization 

I 



and S 3 " (2, 1) is obtained from (p~3|) by i) interchanging the 
label 1 with 2 and ii) replacing the positive sign in the 
denominator with a negative sign. This latter replace- 
ment is given by the fact that, while the second photon 
is emitted, the first photon is absorbed by the atom. Fol- 
lowing the notation used in Refs. [13, HJ, in Eq. (fl~2"j) we 
have furthermore defined 

6*- (1,2) = ^(-l) m ^/ +m ^(2> + l) 1 / 2 (14) 

x ( if L i 3v \ f ju L 2 ji \ 
I —m,j, Ai nij u J \ —m^ M 2 m,j i J 1 

where J "(2,1) is obtained from Eq. (| 14[) by replacing 
L\ O L 2 and Ai O M 2 . 

Equation (TT2]) can be used to evaluate Raylcigh scatter- 
ing by hydrogcnic as well as many-electron systems. For 
the latter case, we must just replace hydrogenic states 
with states of many-electron systems. In this article, re- 
sults are only given for Raylcigh scattering by atomic 
hydrogen. 

C. Evaluation of the transition amplitude for linear 
polarization scenario 

We shall here evaluate the amplitude (fT]) for the case 
in which the incident light is linearly polarized or un- 
polarized and the polarization of the scattered light is 
measured in the linear base. To this end, wc make a de- 
composition of the vector plane wave e l x e tk r by using the 



axis along the z\\k\ direction (as showed in Sec. [TT|) . 
This definition of the quantization axis enables us to 
describe the second photon direction by means of the 
single polar angle 6. Thus the Wigner rotation ma- 
trices simplify as D^ aXa {(p k3 , 6 k2 , 0) = d^ iM {B) and 

Combining Eqs. (JTJ) , ([2]) and (fTTj) . and making use of 
the Wigner-Eckart theorem [36[ , the transition amplitude 
can be written as 



(12) 



(13) 



relation (3|| 

4=4e e_<Ax ^ (is) 

VZ A=±l 

together with Eq. (fTTj) . By combining Eqs. flU, © and 
(|15[) . we get a simple equation which relates the ampli- 
tudes for circular and linear polarization scenarios: 

M\ X i,X2) = \Y, e^ AlXl e 4A « 2 ^ c (A 1; A 2 ) , (16) 

Ai A2 

where for simplicity we defined 

M\ X uX2)=M if {e l xl ,e l X2 ) . (17) 

The amplitude for the linear polarization scenario can be 
thus easily evaluated by using Eqs. (fi"6|) and (fT2)) . 

IV. COMPUTATION 

In this section, we shall discuss how the reduced matrix 
element (fT3|) is calculated in the present work. 

During the last decade, various methods have been 
investigated for calculating the reduced matrix element 
(fT3"|) as well as the whole transition amplitude ((T|) (see 
Refs. [2l| - |23l . [38l l39j ) . In practice, the summation over 
the complete spectrum contained in Eq. (|13[) cannot be 
performed explicitly. Several approaches and approxima- 
tion techniques have thus been proposed in the literature 



_M C (A!,A 2 ) = 2n ]T ^(+*) il - i2+Pl+ ^[ii,i 2 ] 1/2 (A 1 r(A 2 r4| 2 _ A2 (0) 
LiL 2 P1P2 

M 2 

x EC" 1 )" 3 " {2]v I iy/ 2 (^(1,2)5^(1,2) + 6^(2,1)^(2, 1)) , 
where the reduced (second-order) matrix element is given by 



S 3 "(l 2) - V ^'^1" ' tf^Jw )\Pv,3^){Pv,3A a ■ a L 2 ( fc 2,r)||A,Ji> 
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to perform this summation, such as the Coulomb-Green 
function approach, which has been widely used for study- 
ing different decay and scattering processes of atoms and 
ions |UE|]- 

An alternative approach is the finite basis set method 
[42T - |48j . The finite basis set method is based on the sup- 
position that the ion (or atom) is enclosed in a finite 
cavity. Such a restriction leads to a "discretized" contin- 
uum part of the atomic or ionic spectrum, and hence to 
a representation of the Dirac wavefunctions in terms of 
pseudo basis set functions j4^|. The radius of the cav- 
ity R is of course taken large enough to ensure a good 
approximation for the wave functions. 

In the present work, we calculate the reduced matrix 
element (|T5|) by adopting the finite basis set method, us- 
ing B splines and B polynomials as finite basis sets. The 
B splines are one of the most commonly used family of 
piecewise polynomials, since they arc well adapted to nu- 
merical tasks 14311 . The B polynomials, or the Bcrncstein 
polynomials [44j , are a good alternative to the B splines 
since they allow for analytical finite basis-set calculations. 
These are polynomial functions of nth degree that are 
used to obtain the solution of linear and nonlinear dif- 
ferential equations |44|. The details of these basis sets, 
as well as a comparison between them, can be found in 
Ref. j49|. Thus, we restrict ourselves to describe the 
characteristic parameters used in this work. 

The parameters of the B splines basis set are the radius 
of the cavity (i?bs), the number of B splines (ribs) and 
their degree (k). As for the B polynomials, the param- 
eters are the radius of the cavity (i?bp) and the number 
of B polynomials (n\, p ) (the degree of the B polynomials 
is n,b p — 1). The parameters used in both basis sets were 
optimized in order to obtain stability and agreement of 
six digits between the results of both basis sets. The op- 
timal parameters are: i?bs = 60 a.u., ribs = 60, k = 9, 



50 a.u. and ribp = 40. Such set of parameters 



was alre ady obtained for the case of two photon emis- 
sion [ill Hol , and for the angular distribution in Rayleigh 
scattering by atomic hydrogen (30j . 



V. RESULTS AND DISCUSSION 

In this section, we shall present the results for the 
PDADs and the degrees of circular and linear polariza- 
tion of the scattered light which we defined in Sec. IIII Al 
(see Eqs. flU, @, © and ©). Such results have been 
obtained by using the relations presented in Sees. IIII Bl 
and IIII CI and by implementing the computation tech- 
nique presented in Sec. IIVI 



A. Circular polarization scenario 

Fig. [5] displays the angular distribution function (fj|, 
d<7 c /df2, for the four polarization configurations Ai = ±1, 
A2 = ±1 and Ai = ±1, A2 = =pl. The analyzed photon 



energies are 500 eV (a) and 5 keV (6). The shape of 
the angular distribution for E 1 = 500 eV can be easily 
understood owing to the conservation of the angular mo- 
mentum: Since the leading electric-dipole (E1E1) term is 
independent of spin interaction operators within the non- 
rclativistic limit, the photon hclicity must be conserved 
(flipped) for forward (backward) scattering. A similar 
discussion also applies to the higher en ergy case E^ — 
5 keV. However, as discussed elsewhere [301, rnultipolcs 
beyond the E1E1 approximation determine here a strong 
suppression of backwards scattering, thereby suppressing 
helicity-flip scattering events. 

Fig. |3] displays Pc (i.e., the degree of circular polar- 
ization of the scattered light), as a function of the scat- 
tering angle 9, for the two energy values E 1 = 500 eV 
and 5 keV. When compared with Fig. 2, it is evident 
that the degree of circular polarization is less sensitive to 
corrections from higher multipoles than the angular dis- 
tribution. At E 1 = 5 keV, the (weak) effect that higher 
multipoles have on the degree of circular polarization of 
the scattered light, Pc, is to slightly increase it within 
the interval < 9 < 60°, where 90% of the scattering 
events occur pol ]. 

The angular distribution function dtr c /dn can be ob- 
tained from Fig. [2] by taking the arithmetic mean of 
the curves referring to positive (Ai = +1) and negative 
(Ai = —1) helicity of the first photon, for fixed A2. Fol- 
lowing this procedure, the angular distribution takes the 
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FIG. 2: (Color online) Angular distribution of the scattered 
light when the incident light is circularly polarized and the 
scattered light is measured in the circular base (d<7 c /dfi). Re- 
sults are calculated with the account of all photon multipoles 
(solid-black line) and within the electric dipole approximation 
(dashed-red line), for two selected photon energies: a) — 
500 eV; 6) E~,= 5 keV. 
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FIG. 3: (Color online) Degree of circular polarization of the 
scattered light when incident light is circularly polarized (Pc), 
as a function of the scattering angle 0. The curves correspond- 
ing to positive (Ai = +1) and negative (Ai = —1) helicity of 
the incident photon are displayed. Results are calculated with 
the account of all photon multipoles (solid-black line) and 
within the electric dipole approximation (dashed-red line). 



familiar shape ~ 1+cos 2 9 at low energies, independently 
of the helicity of the scattered photon. 

The function Pc (i.e., the degree of circular polariza- 
tion of the scattered light for unpolarizcd incident light) 
vanishes for all angles and energies. This is easily seen 
by applying the definition ^ to the graphs showed in 
Fig. [21 Such a result is somehow expected since, in this 
case, both hydrogen atom and incident light are unpolar- 
ized and therefore there cannot be any preferred direc- 
tion for the circular polarization of the scattered light. If 
there were any, then violation of parity would occur. 



B. Linear polarization scenario 

Fig. U displays the angular distribution function 
dcr'/dO for the four polarization configurations X \ = 
0°(90°), X 2 = 0°(90°) and X i = 0°(90°), X 2 = 90°(0°). 
The analyzed photon energies are 500 eV (a) and 5 keV 
(b). Similarly to the circular polarization scenario, de- 
viations from the electric-dipole approximation are more 
pronounced for higher photon energies. At low energies, 
where the electric-dipole approximation holds, we recover 



Xl 



■e X2 \ , which 



the well-known behavior da (XI1X2) °c |e' 
characterizes the polarization-dependent angular distri- 
bution in the non-rclativistic (Thomson) limit. More 
specifically, at low energies we obtain 



da 1 (0° , 0° ) oc cos^ 9 , da 1 (90° , 90° ) oc (const) , 
do-'(0°,90°) = d<7'(90 o ,0 o ) = 0. 



(18) 



These non-relativistic polarization correlations are the 
same as for two-photon decay in hydrogen atom Q ■ This 
result is not unexpected in view of the similarity of the 
amplitude for the two processes. 

As seen from Fig. [U the angular distribution for 
Rayleigh scattering vanishes at normal angle (9 « 90°) 
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FIG. 4: (Color online) Angular distribution of the scattered 
light when the incident light is linearly polarized and the 
scattered light is measured in the linear base (dcr /dfi). Re- 
sults are calculated with the account of all photon multipoles 
(solid-black line) and within the electric dipole approximation 
(dashed-red line), for two selected photon energies: a) Ey — 
500 eV\ b) E~ - 5 keV. 



if incident light is linearly polarized along xi = 0° direc- 
tion. Moreover, the same figure shows that spin-flip scat- 
tering events (i.e., events for which X i = 0°, X2 — 90° or 
Xi = 90°, X2 = 0°) are strongly suppressed, at any angle 
9 and any energy E y . The linear polarization of the pho- 
ton, which coincides with its spin, is therefore conserved, 
at any photon energy. 

As can be easily verified from Eq. ([7]), suppression of 
spin-flip transitions implies also that the degree of linear 
polarization, Pl, turns out to be simply pa 1, for incident 
polarization along X l — 0°, 90° direction, and therefore 
is not displayed. 

The results for the function da 1 /dfi can be obtained 
from Fig. 01 by taking the arithmetic mean of the curves 
referring to X i = 0° and xi = 90°, for fixed X2- 
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FIG. 5: (Color online) Degree of linear polarization of the 
scattered light when incident light is unpolarized (Pl), as a 
function of the scattering angle 9. Results are calculated with 
the account of all photon multipoles (solid-black line) and 
within the electric dipole approximation (dashed-red line). 

Fig. [5] displays Pl (i.e, the degree of linear polariza- 
tion of the scattered light for unpolarized incident light) 
as a function of the scattering angle B, for the two pho- 
ton energies 500 eV and 5 keV. We read from the figure 
that light scattered at normal angle is fully linearly po- 
larized, while light scattered forwards or backwards is 
unpolarized, at any photon energy. This polarization 
feature is also known to characterize the scattering of 
light by small molecules, when the photon wavelength is 
small compared with the molecular radius. It is due to 
this property that the sunlight scattered by surfaces at 
normal angles is always linearly polarized along the axis 
which is orthogonal to the scattering plane Q ■ 

VI. SUMMARY AND PERSPECTIVES 

We studied polarization properties of light in Raylcigh 
scattering by unpolarized atomic hydrogen, based on 
second-order perturbation theory and Dirac relativistic 
equation. We introduced two experimental scenarios in 
which the light is circularly and linearly polarized, re- 
spectively. For each of these scenarios, we analyzed the 
polarization-dependent angular distribution and the de- 
grees of circular and linear polarization of the scattered 
light. To this end, we first decomposed the Rayleigh 
transition amplitude in terms of spherical tensors (an- 
gular part) and reduced amplitudes (radial part). We 
then calculated these latter for scattering by hydrogen 
atom by means of the finite basis-set method based on 



the relativistic Dirac equation. The polarization depen- 
dent angular distribution and the degrees of circular and 
linear polarization of the scattered light were plotted for 
the photon energy range 0.5 to 10 keV. 

We found that, for circularly polarized incident light, 
hclicity-flip scattering events are allowed at low energies 
(Ey < 500 eV) and are suppressed at high energies (E y > 
5 keV), due to the suppression of backward scattering. 
Thus, the helicity of the incident photon is not conserved 
at low photon energies while is conserved at high photon 
energies. 

For linearly polarized incident light, it was showed that 
spin-flip scattering events (i.e., events for which the linear 
polarization angle of incident and scattered light differ by 
a normal angle) are strongly suppressed, at any scattering 
angle and photon energy. Thus, the linear polarization of 
the incident photon is approximately conserved and con- 
veyed to the scattered photon, in the whole investigated 
photon energy range. 

Finally, light scattered at normal angles was found to 
be fully linearly polarized, at any photon energy value. 

With advances in polarization sensitive detectors in the 
x-ray and 7-ray region, theoretical studies on polarization 
of photons emitted or scattered by atomic systems have 
lately become important 0, [l5|, [l8| ■ Further studies to 
investigate polarization properties of light scattered by 
solid targets are underway. 



VII. ACKNOWLEDGMENT 

L. S. and F. F. acknowledge financial support by the 
Research Council for Natural Sciences and Engineer- 
ing of the Academy of Finland. P. A. and S. T. ac- 
knowledges the support of German Research Founda- 
tion (DFG) within the Emmy Noether program under 
Contract No. TA 740 1-1. S. F. acknowledges sup- 
port by the FiDiPro program of the Finnish Academy. 
J. P. S. and P. A. acknowledge the support by FCT 
- Fundacao para a Ciencia e a Tecnologia (Portugal), 
through the Projects No. PEstOE/FIS/UI0303/2011 
and PTDC/FIS/117606/2010, financed by the European 
Community Fund FEDER through the COMPETE - 
Competitiveness Factors Operational Programme. F.F. 
acknowledges financial support by Fundagao de Amparo 
a Pesquisa do estado de Minas Gerais (FAPEMIG). 



[1] L. Sharma, A. Surzhykov, M. K. Inal, and S. Fritzsche, 

Phys. Rev. 81, 023419 (2010). 
[2] A. Bondarevskaya, L. Labzowsky, A. Prozorov, G. Plu- 

nien, D. Liesen, and F. Bosch, J. Phys. B 43, 245001 

(2010). 

[3] F. Fratini and A. Surzhykov, Hyp. Int. 199, 85 (2011). 
[4] F. Fratini, M. C. Tichy, Th. Jahrsetz, A. Buchleitner, 



S. Fritzsche, and A. Surzhykov, Phys. Rev. A 83, 032506 
(2011). 

[5] S. Tashenov, T. Back, R. Barday, B. Cederwall, J. En- 
ders, A. Khaplanov, Yu. Poltoratska, K.-U. Schassburger, 
and A. Surzhykov, Phys. Rev. Lett. 107, 173201 (2011). 

[6] S. Tashenov, et al, Phys. Rev. Lett. 97, 223202 (2006). 

[7] C. F. Bohren and A. B. Fraser, Phys. Teach. 23, 267 



8 



(1985). 

[8] P. P. Kane, L. Kissel, R. H. Pratt, and S. C. Roy, Phys. [29 
Rep. 140, 75 (1986). 

A. Boke, Radiat. Phys. Chem. 80, 609 (2011). [30 
N. S. Kampel, A. Griesmaier, M. P. Hornbak Steenstrup, 

F. Kaminski, E. S. Polzik, and J. H. Miiller, Phys. Rev. [31 
Lett. 108, 090401 (2012). 

D. G. Norris, A. D. Cimmarusti, L. A. Orozco, [32 
P. Barberis-Blostein, and H. J. Carmichael, Phys. Rev. 
A 85, 021804 (2012). [33 
L. V. Kulik, K. Ovchinnikov, A. S. Zhuravlev, V. E. Bisti, 
I. V. Kukushkin, S. Schmult, and W. Dietsche, Phys. [34 
Rev. B 85, 113403 (2012). 

V. Bellani, F. Rossella, M. Amado, E. Diez, K. Kowalik, 

G. Biasiol, and L. Sorba, Phys. Rev. B 83, 193307 (2011). 
J. Crassous, A. Anion, and J. Crassous, Phys. Rev. A [35 
85, 023806 (2012). 

S. Tashenov, A. Khaplanov, B. Cederwall, and K.-U. [36 
Schassburger, Nucl. Instrum. Methods Phys. Res. Sect. 
A 600, 599 (2009). [37 
C. Spezzani, E. Allaria, M. Coreno, B. Diviacco, E. Fer- 
rari, G. Geloni, E. Karantzoulis, B. Mahieu, M. Vento, [38 
and G. De Ninno, Phys. Rev. Lett. 107, 084801 (2011). 
F. Smend, D. Schaupp, H. Czerwinski, M. Schumacher, [39 
A. H. Millhouse, and L. Kissel, Phys. Rev. A 36, 5189 
(1987). 

A. Bondarev, GSI Summer Student Program 2009: In- [40 
ternal Report, Rayleigh scattering in X-ray polarimetry; 
A. Gumberidze and M. Schwemlein, private communica- [41 
tion. 

M. Gavrila, Phys. Rev. 163, 147 (1967). [42 

H. R. Sadeghpour and A. Dalgarno, J. Phys. B 25, 4801 
(1992). [43 
H. M. T. Nganso and M. G. K. Njock, J. Phys. B 40, 807 
(2007). [44 
A. Costescu, K. Karim, M. Moldovan, S. Spanulescu, 
and C. Stoica, J. Phys. B 44, 0452040 (2011). [45 
A. Costescu, P. M. Bergstrom, Jr., C. Dinu, and 
R. H. Pratt, Phys. Rev. A 50, 1390 (1994). [46 
A. Costescu, S. Spanulescu, and C. Stoica, J. Phys. B 
40, 2995 (2007). [47 
J. P. J. Carney, R. H. Pratt, N. L. Manakov, and A. V. 
Meremianin, Phys. Rev. A 61, 042704 (2000). [ 
S. C. Roy, B. Sarkar, L. D. Kissel, and R. H. Pratt, Phys. 
Rev. A 34, 1178 (1986). [49 
N. L. Manakov, A. V. Meremianin, A. Maquet, and J. P. 
J. Carney, J. Phys. B 33, 4425 (2000). [50 
N. L. Manakov, A. A. Nekipelov, and A. G. Fainshtein, 
Yad. Fiz. 45, 1091 (1987) [Sov. J. Nucl. Phys. 45, 677 



(1987)]. 

N. L. Manakov, A. V. Meremianin, J. P. J. Carney, and 
R. H. Pratt, Phys. Rev. A 61, 032711 (2000). 
L. Safari, P. Amaro, S. Fritzsche, J. P. Santos, and F. 
Fratini, Phys. Rev. A 85, 043406 (2012). 
F. Fratini, One- and Two-photon decays in Atoms and 
Ions (LAP, Saarbriicken, Germany, 2011). 
A. I. Akhiezer and V. B Berestetskii, Quantum Electro- 
dynamics (Wiley, New York, 1965). 

F. Fratini and A. G. Hayrapetyan, Phys. Scr. 84, 035008 
(2011). 

V. V. Balashov, A. N. Grum-Grzhimailo and 
N. M. Kabachnik, Polarization and Correlation Phe- 
nomena in Atomic Collisions (Kluwer Academic Plenum 
Publishers, New York 2000). 

M. E. Rose, Elementary Theory of Angular Momentum 
(John Wiley, New York, 1957). 

J. J. Sakurai, Modern Quantum Mechanics (Addison- 
Wesley, 1994). 

S. P. Goldman and G. W. F. Drake, Phys. Rev. A 24, 
183 (1981). 

F. Fratini, S. Trotsenko, S. Tashenov, Th. Stohlker, and 
A. Surzhykov, Phys. Rev. A 83, 052505 (2011). 
A. Surzhykov, A. Volotka, F. Fratini, J. P. Santos, P. 
Indelicato, G. Plunien, Th. Stohlker, and S. Fritzsche, 
Phys. Rev. A 81, 042510 (2010). 

A. Surzhykov, P. Koval, and S. Fritzsche, Phys. Rev. A 
71, 022509 (2005). 

A. Maquet, V. Veniard, and T. A. Marian, J. Phys. B 
31, 3743 (1998). 

J. Sapirstein and W. R. Johnson, J. Phys. B 29, 5213 
(1996). 

W. R. Johnson, S. A. Blundell, and J. Sapirstein, Phys. 
Rev. A 37, 307 (1988). 

D. D. Bhatta and M. I. Bhatti, Appl. Math. Comput. 
174, 1255 (2006). 

O. Zatsarinny and K. Bartschat, Phys. Rev. A 77, 062701 

(2008) . 

C. Froese Fischer and O. Zatsarinny, CPC 180, 879 

(2009) . 

A. Surzhykov, P. Indelicato, J. P. Santos, P. Amaro, and 
S. Fritzsche, Phys. Rev. A 84, 022511 (2011). 
J. P. Santos, F. Parente, and P. Indelicato, Eur. Phys. J. 
D 3, 43 (1998). 

P. Amaro, A. Surzhykov, F. Parente, P. Indelicato, and 
J. P. Santos, J. Phys. A 44, 245302 (2011). 
P. Amaro, J. P. Santos, F. Parente, A. Surzhykov, and 
P. Indelicato, Phys. Rev. A 79, 062504 (2009). 



